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ABSTRACT 


Site  diversity  is  usually  proposed  as  a  means  for  improving 
earth/satellite  path  availability  at  EHF,  because  intense  rain  cells  are 
generally  small  in  diameter.  However,  there  are  many  applications,  such  as 
point-to-point,  high-data-rate  communications,  where  site  diversity  is  either 
very  inconvenient  or  impossible.  In  such  cases,  provided  the  user  can  accept 
a  short  delay,  Incoming  data  can  be  stored  and  retransmitted  at  a  higher  rate 
once  the  period  of  high  attenuation  has  passed.  By  using  statistical  long¬ 
term  attenuation  and  fade-duration  probability  distributions  gathered  at 
several  frequencies  and  geographic  locations,  an  empirical  model  for 
determining  the  improvement  in  effective  availability  created  by  a  given 
amount  of  storage  time  is  developed.  The  model  appears  to  be  applicable  to 
all  geographical  areas  where  standard  climatological  data  is  gathered. 
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I.  INTRODUCTION 


1.1  Rationale 

There  are  many  applications  where  site  diversity  is  either  difficult  or 
Impossible  for  small  EHF  terminals.  The  difficulty  may  be  due  to  factors  such 
as  cost,  tactics  (battlefield  communications  may  preclude  the  interlinking  of 
terminals  by  other  communications  capabilities)  or  geographic  constraints  (it 
may  not  be  possible  to  acquire  real  estate  the  appropriate  distance  away  from 
a  given  terminal  location).  In  such  cases,  it  may  be  useful  to  store  data 
during  periods  of  high  rain-induced  path  attenuation  and  transmit  at  a 
somewhat  higher  rate  during  the  next  period  of  low  attenuation. 

The  availability  of  an  earth-satellite  path,  PA,  is  defined  as  one  minus 
the  fraction  of  the  time  over  an  average  year  that  the  system  path  margin  is 
insufficient  to  maintain  the  desired  link  probability  of  error.  In  this 
report  only  rain  is  assumed  to  affect  the  system  path  margin;  other  factors, 
such  as  component  failure,  atmospheric  absorption  repair  time,  suntranslt,  and 
antenna  pointing  errors  are  not  considered.  If  the  margin  is  M  (dB),  the 
probability  of  rain  is  PraIN*  and  the  c“®“l*tive  probability  density 
distribution  of  attenuation,  A,  when  it  is  raining,  is  P(A>AQ| given  rain),  the 
system  availability,  for  AQ«M,  becomes 

P  -  1  -  P(A  >  A  )  -  1  -  P_.TU  •  P( A  >  Mlgiven  rain)  (1.1.1) 

A  O  KAJ.N 

There  are  numerous  studies  which  attempt  to  link  P(A>AQ)  with  known  rain- 
rate  distributions  and  climatological  factors  [see,  for  example, 
Engelbrecht,  1979;  Lin,  1973;  Bussey,  1950;  Hogg  and  Chu,  1975;  Rogers,  1976; 
Crane  and  Blood,  1979;  Kaul,  Wallace  and  Klnal,  1980;  Ippolito,  1978;  and 
Crane,  1971].  This  study  does  not  attempt  such  linkage,  except  where 
absolutely  necessary*,  but,  instead,  looks  to  see  if  any  "universal”  behaviors 
can  be  found. 


* 


or  where  the  author's  interest  was  piqued. 
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In  Section  2.1  the  basic  model  elements  for  the  attenuation  and  fade** 
duration  distributions  are  presented.  A  model  for  the  probability  of  rain  is 
developed  in  Section  2.2.  All  of  the  necessary  climatologic  data  for  the 
locations  and  time  spans  for  which  good  attenuation  data  is  available  are 
given  in  Section  2.3.  Section  2.4  contains  a  rationale  for  Ignoring 
elevation-angle  dependencies  and  Section  2.5  examines  the  Fresnel-zone- 
influenced  time  structure  of  rain  fading.  The  measured  data  for  several 
different  sites  are  detailed  in  Section  2.6  and  the  basic  inputs  to  a  model 
for  storage-time  effects  are  developed  for  use  in  the  storage-time  model 
presented  in  Section  2.7.  Finally,  an  example  of  the  use  of  the  storage-time 
model  is  given  in  Section  2.8. 

1.2  Conclusions 

The  basic  conclusions  are  given  in  Section  III  of  this  report. 
Essentially,  the  effective  margin  M2  (dB)  with  storage  time  of  T  (minutes),  in 
terms  of  the  margin  Mj^  (dB)  without  storage,  becomes 

T  3.7D  /D  ]0.5 

M2/M1  -  [l  +  0.2  T'e  1  KJ  (1.2.1) 

where  ftp  is  the  number  of  thunderstorm  days  per  year  and  Dg  is  the  number  of 
days  per  year  of  rain  (as  defined  in  Section  2.3).  This  model  appears  to  hold 
for  all  climatologic  regions  of  the  globe.  For  a  climate  similar  to  that  of 
eastern  Massachusetts,  where  Dp-19  and  DR-117,  a  storage  time  of  10  minutes 
will  provide  an  effective  increase  in  the  path- loss  margin  of  2.16  times. 
That  is,  the  system  availability  at,  say,  10  dB  is  improved  to  that  that  would 
have  existed  with  a  margin  of  21.6  dB. 


II.  ANALYSIS 


2.1  Model  Elements  -  The  Lognormal  Distribution 

If  attenuation  data  Is  taken  over  a  long  period  of  time,  generally  in 
excess  of  1  year,  the  distribution  of  attenuation  during  the  time  it  is 
raining  will  be  found  to  be  lognormal.  This  lognormal  behavior  was  recognized 
by  Lin  [1973]  and  recently  reviewed  by  Engelbrecht  [1979].  In  addition,  in 
most  cases  the  distribution  of  the  length  of  time  that  the  attenuation  exceeds 
a  given  threshold  is  also  lognormally  distributed  [Lin,  1973].  Specifically, 
the  attenuation  in  dB  will  be  distributed  as 


where 


(2.1.1) 


P(A>Aq)  is  the  probability  that  the  path  attenuation  A  (dB)  exceeds 
AQ(dB), 

PraIN  1s  the  Pr°bability  that  it  is  raining, 

An  is  median  value  of  the  path  attenuation  (dB)  during  the  time 
It  is  raining,  and 

S^  is  the  standard  deviation  of  InA  during  the  time  it  is 
raining. 


For  a  lognormal  distribution  the  mean  is  related  to  the  median  as 


A  =* 


A  e 
ra 


,/2 


dB 


(2.1.2) 


In  much  the  same  fashion,  the  cumulative  distribution  of  fade  durations 
beyond  a  given  threshold,  AQ(dB),  becomes 


P(T 


>  V 


(2.1.3) 
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where 


P(x>xQ)  is  the  probability  that  the  fade  duration  t  (minutes)  exceeds 
xQ  (minutes), 

t  is  the  median  fade  duration  in  minutes,  and 
m 

ST  is  the  standard  deviation  of  lnx. 


The  median  value  of  x  is 


i 


St/2 


x  e 
m 


and  it  should  be  realized  that  x,  x„,  and  S,.  are  functions  of  the  chosen 

'  m*  t 

attenuation  threshold,  AQ. 

Now,  the  probability  that  the  rain  attenuation  exceeds  a  given  threshold 
AQ  for  a  time  in  excess  of  t0(A0)  becomes  the  product  of  3  terms 


RAIN 


2  e*fc 


InA  -InA  ,  flnx  (A  )-lnx  (A  ) 

o  m  1  ,  oo  mo,-, 

-  •  j  erfc  - — -  (2.1. 

/2  S,  z  I  /2  S  (A  ) 

(.  A  "*  L  T  o  J 


5) 


It  is  hoped  that  most  of  the  non-time-stationary  behavior  of  rain-induced 


attenuation  can  be  represented  by  the  first  term,  P, 


RAIN* 


and  that  the  next  two 


terms  can  be  considered  stationary.  If  so,  a  system  with  no  data  storage 
capability  and  a  margin  M^(dB)  will  have  a  probability  of  exceeding  during 
a  rainstorm 


P(M  >  Mx) 


(2.1.6) 


If  a  storage  time  of  xQ  is  available  then  a  new,  equivalent  margin,  h^dB), 
can  be  defined  from  (2.1.5)  and  (2.1.6)  as 


1 

T 


j  erfc 


lnx  (M, )-lnx  (M,  ) 
o  1 _ m  1 

/2  S 

x 

(2.1.7) 
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A  simple  relation  is  sought  for  the  increase  in  effective  margin, 

The  "increase"  in  margin  is  really  a  reflection  of  the  increase  in 
availability  gained  by  storage  and  discussed  by  Engelbrecht  [1979].  Generally 
M2^M1  be  a  function  of  A^,  S^,  Tm,  ST,  and  .  However,  as  will 
be  seen,  the  explicit  dependence  on  all  but  t0(mi)  is  slight. 

Because  the  path-attenuation  statistics  are  complex  functions  of  rainfall- 
rate  statistics  and  such  rainfall  statistics  are  functions  of  the  fraction  of 
the  total  rainfall  contained  in  thunderstorms  (convective  rain),  it  should  be 
expected  that  will  have  some  functional  dependence  on  the  probability  of 
occurance  of  thunderstorms.  Such  a  potential  dependence  is  strengthened  by 
the  correspondences  found  in  Section  2.3.  In  Section  2.7  will  be  found 
to  be  a  function  of  the  ratio  of  days/year  with  thunderstorms  to  days/year 
with  rain. 

2.2  On  the  Determination  of  the  Probability  of  Rain 

Some  of  the  attenuation  data  available  in  the  literature  is  presented  as  a 
percent  of  the  year  a  given  attenuation  is  exceeded  and  some  is  presented  as 
percent  of  the  raining  time  a  given  attenuation  is  exceeded.  In  order  to 
convert  from  one  format  to  the  other,  the  probability  that  it  will  rain  must 
be  known.  Such  a  probability  of  rain  also  places  a  lower  limit  on  rain- 
induced  availability;  that  is,  a  PraIN”®*^  implies  an  availability  always 
greater  than  0.95. 

The  value  of  represents  an  instantaneous  measure  of  the  percentage 
of  time  it  rains.  Unfortunately,  the  available  climatological  data  seldom 
Includes  measures  of  instantaneous  rain.  The  usual  statistics  given  include 
such  measures  as  yearly  totals  of  precipitation,  thunderstorm  days,  cloudy 
days,  total  snowfall  and  number  of  days  when  precipitation  exceeded  0.25  mm 
(0.01  inch). 

Correlations  between  percentage  of  time  a  given  rain  rate  is  exceeded  and 
the  ratio  of  annual  precipitation  to  number  of  days  with  measurable 
precipitation  (0.25  mm)  have  been  observed  [see  Handbook  of  Geophysics, 
1961].  However,  the  rain  rates  are  clock-hour  rates  and  not  instantaneous 
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values  (for  example,  minute-by-minute  or  second-by-second)  and  again  the 
percentages  are  of  a  total  year  not  the  raining  time.  Other  means  must  be 
found  to  determine 

Bussey  [1950],  in  one  of  the  first  comprehensive  treatments  of  the  effect 
of  rain  on  EHF  path  attenuation,  determined  relationships  between 
instantaneous,  1 / 2— ,  1-,  and  2-hour-interval  measures  of  rain  rate  for 
Washington,  DC.  His  "instantaneous"  measures  were  1-minute  rates.  The  longer 
the  interval  over  which  rainfall  is  determined,  the  higher  the  percent  of  the 
year  it  would  be  effectively  raining.  That  is,  if  a  minute-by-minute 

determination  of  rain/no-rain  is  kept,  the  percentage  of  raining  minutes  out 
of  the  total  minutes  in  a  year  would  be  a  smaller  figure  than  would  be  yielded 
by  a  similar  determination  of  day-by-day  rain.  The  reason  for  this  difference 
is  obvious:  most  days  during  which  rain  falls  are  not  "complete”  and  many 

will  have  only  a  small  fraction  of  the  day  with  rainfall;  however,  most  rain 

* 

minutes  are  not  so  fractionated  . 

The  long-term  precipitation  data  for  Washington,  DC  (taken  as  the  average 
of  National-  and  Dulles-Airport  measurements)  is  presented  in  Table  2.2.1. 
Although  the  water  content  of  snow  varies  widely  from  values  probably  less 
than  5%  in  northern  Alaska  in  January  to  as  great  as  20%  in  the  "slushy" 

spring  snowfalls  common  along  the  east  coast  of  the  US,  an  arbitrary  figure  of 
10%  was  chosen  to  remove  snowfall  from  the  total  precipitation  data  (since 
snow  attenuation  is  much  smaller  than  that  produced  by  rain  even  at  the  same 

"water"  rates).  Similarly,  in  order  to  estimate  the  number  of  days  of 

measurable  rainfall,  the  number  of  days  with  snow  greater  than  25.4  mm  were 
subtracted  from  days  of  measurable  precipitation  (perhaps  days  of  snow  greater 
than  2.5  mm  would  have  been  more  appropriate,  but  such  a  measure  is  not 
available  in  the  National  Summarys  of  Climatological  Data). 

If  Bussey's  {1950]  relationships  are  examined  at  the  very  lowest  rain 
rates  (and  extrapolated  somewhat  for  the  1/2-,  1-,  and  2-hour  interval  values) 

*A  more  recent  analysis  of  the  conversion  from  one  averaging  time  to  another 
was  presented  by  Lin  [1976a]. 


and  If  the  Instantaneous  hours/year  are  set  at  500  (for  a  1-minute  rain  rate 
of  0.13  mm/hour  or  a  total  of  2xlO“3  mm  in  1-rainute),  the  ratio,  r,  of  1/2-, 
1-,  2-  and  24-hour  values  to  the  instantaneous  value  can  be  determined.  Using 
the  107-day  figure  from  Table  2.2.1,  the  24-hour  value  becomes  107/365.25  or 
0.293.  The  instantaneous  value  is  0.057  (500  hours  in  8760).  The  24-hour 
value  and  the  intermediate  values  of  r  are  plotted  in  Fig.  2.2.1. 

The  simple  power-law  relation  for  r-1  is  probably  representative  of  many 
different  climates.  As  a  result,  if  the  fraction  of  the  year  given  by 
0R/365.25  is  divided  by  the  factor  of  5.14  from  Fig.  2.2.1,  a  reasonable  value 
for  PRA^N  can  be  obtained.  That  is,  if  the  number  of  days  of  measurable 
(>0.25  mm)  precipitation  are  reduced  by  the  number  of  days  of  snow  greater 
than  25.4  mm  and  the  figure  that  results  is  multiplied  by  the  reciprocal  of 
365.25x5.14,  the  resulting  probability  of  rain  can  be  written  as 


°p-  °S 


RAIN  (365.25X5.14) 


5.3  x  10-4  D„ 


(2.2.1) 


It  is  expected  that  such  a  derived  value  of  PRAIN  should  be  reasonably  valid 
for  climates  as  widely  varied  as  Yuma,  AZ,  with  D^lb  and  pra.IN=0*00®3  to 
Hilo,  HI,  with  DR=282  and  PraIN*®'1 3 ’ 

Lin  [1975]  nas  defined  PRAIN  from  1-hour  integration-time  data  supplied  by 
the  Weather  Bureau.  His  values  of  Prain  are  1-minute  averages  set  equal  to 
1/2  of  the  1-hour  rain  rates  greater  than  0.25  mm/hr.  The  values  he  records 
for  different  locations  compare  favorably  with  those  using  the  DR  approach 
developed  here.  In  the  same  paper,  Lin  [1975]  also  presents  an  empirical 
formula  which  increases  the  value  of  Pradj  obtained  from  point  rain  data  to 
account  for  the  Increased  probability  of  occurance  of  rain  along  a  path  over 
that  expected  at  a  point.  His  formula  is  based  on  Florida  thunderstorm  data 
for  terrestrial  paths  and  shows  quite  substantial  increases  for  moderate 
length  paths  (~10  km  would  Increase  a  point  value  of  Prain*0*060  t0  0.083). 
For  an  area  like  Florida  with  a  high  rate  of  scattered  thunderstorms  such  an 
increase  might  well  be  appropriate.  However,  because  of  the  natural  movement 
of  thunderstorms,  the  increase  should  be  azimuthly  sensitive,  Lin’s  is  not; 


TABLE  2.2.1 

LONG-TERM,  NOMINAL  PRECIPITATION  DATA  FOR  WASHINGTON,  D.C.* 


WT,  Total  Precipitation  In  ram  1004 

Wg,  Total  Snowfall  In  mm  493 

WR,  Total  Rainfall  In  ram 

(taken  as  WT  -  0.1  Wg)  955 

Dp,  Number  of  Days  Precipitation  >  0.25  mm  112.5 

Dg,  Number  of  Days  Snowfall  >  25.4  mm  5.5 


DR,  Number  of  Days  of  Measurable  Rain 
(taken  as  Dp  -  Dg)  107 

DT,  Number  of  Thunderstorm  Days  28 

*Data  are  averages  of  National  and  Dulles  Airport  measurements. 
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moreover,  the  elevated  paths  for  earth/satellite  communications  should  behave 
differently  than  terrestrial  paths  and  other  climate  areas  with  a  smaller 
proportion  of  rain  in  thunderstorms  should  show  less  of  an  increase. 
Consequently,  such  a  refinement  will  not  be  attempted  here  and  the  value  of 
PRAIN  ^or  a  P°^nt»  derived  from  D^,  will  be  used  in  the  following  analyses. 

Finally,  although  the  values  of  P^^  and  relative  thunderstorm  index, 
Dt/Dr,  where  DT  is  the  number  of  yearly  thunderstorm  days,  are  long-term 
average  values,  there  can  be  substantial  short-term  variability  in  both 
measures.  If,  for  some  reason,  knowledge  of  the  seasonal  variability  of  PRAIN 
and  D^/Dr  is  important,  a  monthly  variation  can  be  computed.  For  a  monthly 
value  of  ca^^ed  PRAxn  M*  (2.2.1)  can  be  changed  to 

D. 


'PM  ~  °SM 


RAINjM  (days  in  that  month)(5. 14) 


(2.2.2) 


where  Dp;i{  is  the  number  of  days  during  that  month  where  precipitation  greater 
than  0.25  ram  occurred  and  DgM  is  the  number  of  days  during  that  month  where 


snowfall  greater  than  25.4  mm  occurred. 
(2.2.2)  becomes 


If  Drm  is  defined  as  DpM_DSM,  then 


RM 


RAIN,M  (days  in  that  month)(5.14) 


(2.2.3) 


In  similar  fashion,  the  monthly  thunderstorm-day  index,  DTM»  yields  the 
relative  value 

As  an  illustration  of  the  potential  climatological  variations  in  Pr^xn  M 
and  Dtm/Drm  over  the  year,  average  monthly  figures  for  the  period  1941  to  1970 
were  obtained  from  Ruffner  [1978].  The  monthly  variations  for  two  North-East- 
Coast-US  locations  are  plotted  in  Fig.  2.2.2.  Two  representative  Gulf-Coast 
locations  are  plotted  in  Fig.  2.2.3  and  two  North-West-Coast  locations  are 
given  in  Fig.  2.2.4.  A  sample  of  a  Central-Plains  state  and  a  Rocky-Mountain 
state  is  given  in  Fig.  2.2.5  and  two  climate  extremes  for  the  US,  Yuma,  AZ  and 
Hilo,  HI,  are  illustrated  in  Fig.  2.2.6. 
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-Plains  location  and  a  Rocky-Mountain  location 


Several  Interesting  observations  can  be  made.  The  yearly  variability  of 
PRAIN  M  on  tlie  Nortll~East“Coa3t  an(*  the  desert  (Yuma,  AZ)  and  Pacific-Island 
(Hilo,  HI)  locations  is  small.  "April  showers"  may  "bring  May  flowers"  only 
in  Hew  England.  Gulf-Coast  locations  have  a  peak  in  P^j^  ^  in  the  summer 
months  (and  a  second,  minor  peak  in  late  winter)  with  the  summer  peak  the 
result  of  much  increased  thunderstorm  activity.  On  the  other  hand,  the  North- 
West-Coast  locations  have  their  "wet  seasons"  in  late  fall  and  winter  and  the 
summer  months  are  dry. 

The  values  of  that  result  from  an  average  over  the  year  of  PRAIN  M 

for  the  locations  chosen  will  be  slightly  different  from  those  used  later  in 
Section  2.3.  Similar  variations  in  Dtj,/DR  will  be  found.  The  reasons  for 
these  differences  lie  in  the  different  data  bases  used  (1941-1970  versus  a 
much  longer-term  average).  Also,  because  thunderstorm  days  are  defined  as 
days  in  which  thunder  is  heard,  on  a  monthly  basis  dtm^DRM  can  be  greater  than 
unity.  This  behavior  will  be  most  pronounced  in  such  locations  as 
New  Orleans,  LA,  Tampa,  FL,  Denver,  CO,  and  Yuma,  AZ,  where  most  of  the  summer 
rainfall  occurs  as  thunderstorms  and  at  any  given  point  a  thunderstorm  may 
well  be  heard  but  no  rain  falls  from  that  storm  on  that  day  at  that  particular 
location. 

2 . 3  Climatological  Data 

The  best  and  most  complete  long-term  data  available  on  rain-induced  earth- 
space  path  attenuation  at  EHF  comes  from  sites  in  Waltham,  MA  [Nackoney,  1979, 
and  Tang,  1979],  Austin,  TX  [Vogel,  1979],  Crawford  Hill,  NJ  [Wilson,  1969  and 
Lin,  1973],  Tampa,  FL  [Tang,  1979,  and  Tang,  1980],  Slough,  ENG  [Davis,  1971, 
and  Davis,  1973],  and  Rosman,  NC  [Ippolito,  1975  and  Crane,  1980].  In  order 
to  interpret  the  data,  the  climatological  data  for  each  area,  both  nominal  and 
for  the  period  the  data  were  taken,  are  needed. 

The  data  for  the  US  sites  can  be  obtained  from  the  Climatological  Data 
reports  of  the  National  Oceanic  and  Atmospheric  Administration  (sometimes 
called  LCD's  for  Logical  Climatological  Data).  The  climatological  data  for 
Slough,  ENG,  were  obtained  from  several  sources,  including  the  USAF 
Environmental  Technical  Applications  Center  Worldwide  Airfield  Climatic  Data 


(April  1971]  and  the  Handbook  of  Geophysics  [1961],  and  monthly  data  were 

summed  to  obtain  values  which  cover  periods  different  from  one  calendar 
year . 

Nominal  values  for  total  yearly  rainfall,  days  with  rainfall  In  excess  of 
0.25  mm,  thunderstorm  days  and  probability  of  rain  computed  using  (2.2.1)  are 
given  in  Table  2.3.1.  The  same  parameters  for  the  observational  periods  are 
given  In  Table  2.3.2.  It  is  of  Interest  to  note  that  the  Crawford  Hill,  NJ, 
PrAIN  va^ue  of  0.056  for  1968,  which  is  obtained  from  Newark  Airport 

measurements,  compares  favorably  with  the  value  of  0.0523  measured  by  the 
experimenters  and  given  by  Lin  [1972].  In  addition,  Lin  [1972]  gives  the 
value  of  Prajn  ^or  Slough,  ENG,  as  0.085  which  compares  favorably  with  the 
nominal  value  of  0.089  and  the  period  value  of  0.092  as  obtained  by  this 

simple  modification  of  daily  climatological  data. 

It  is  interesting  to  examine  the  instantaneous  rain-rate  curves  redrawn 
taking  P^iN  (derived  from  (2.2.1)  and  climatological  data)  into  account. 
Four  representative  (Middle-East-Coast,  Gulf-Coast,  Central-Plains,  and  North- 
West-Coast)  US  areas  are  chosen  to  exhibit  medium,  very-high,  high  and  low 
thunderstorm  activities,  respectively.  Data  on  rain  rates  from 

Washington,  DC,  New  Orleans,  LA,  Oklahoma  City,  OK,  and  Portland,  OR,  were 
obtained  from  the  Handbook  of  Geophysics  [1961]  as  percent-of-year  values  and 
changed  to  cumulative  probabilities  during  periods  of  rain  using  PRAjN  from 
(2.2.1)  and  nominal  climatological  data.  The  results  are  plotted  in 
Fig .  2.3.1. 

Note  that  the  low-rate  portions  of  the  curves  are  essentially  lognormal 
and  the  high-rate  portions  represent  a  second  distribution  largely  due  to 
thunderstorm  activity*.  Even  Portland,  OR,  with  D^=7,  shows  some  bimodal 
character.  The  fall-off  at  very  high  rates  can  be  due  to  lack  of  sufficient 
statistics  and/or  the  1-minute  averaging  times.  Lin  (1976b)  has  proposed  a 
means  of  estimating  the  extreme-value  statistics.  However,  if  the  averaging 

*This  difference  was  also  noted  by  Hogg  [1969]. 


TABLE  2.3.1 

NOMINAL  CLIMATOLOGIC  DATA  FOR  DATA-GATHERING  SITES 


NOMINAL  VALUES 


LOCATION 

WR,  MMb 

dr,  days' 

dt,  DAYS 

PRAIN 

Waltham,  MA1 2 

971 

117 

19 

0.062 

2 

Rosraan,  NC 

1149 

119 

46 

0.063 

Austin,  TX 

823 

82 

41 

0.044 

Crawford  Hill,  NJ3 4 

981 

115 

25 

0.061 

Tampa ,  FL 

1254 

107 

88 

0.057 

A 

Washington,  DC 

955 

107 

28 

0.057 

Slough,  ENG5 6 7 8 

605 

167 

15 

0.089 

1.  Boston  (Logan  Airport),  MA 

2.  Average  of  Greenville,  SC  and  Asheville,  NC 

3.  Newark  (Newark  Airport),  NJ 

4.  Average  of  National  and  Dulles  Airport  data 

5.  London  (Heathrow  Airport),  ENG 

6.  Total  precipitation  minus  one-tenth  of  total  snow,  sleet,  hail,  etc. 

7.  Number  of  days  precipitation  is  greater  than  0.25  mm  minus  number  of  days 
snow,  sleet,  is  greater  than  25.4  mm 

8.  Taken  as  5.07xl0-^  times  DR 


f 


TABLE  2.3.2 

PERIOD-RELATED  CLIMATOLOGICAL  DATA 


PERIOD  DATA 


LOCATION 

PERIOD 

WR,  MM 

°R‘ 

DAYS 

DT,  DAYS 

PRAIN 

Waltham,  MA1 2 

2/77-1/78 

1065 

(971) 

124 

(117) 

22  (19) 

0.066 

(0.062) 

Waltham,  MA1 

2/78-1/79 

855 

(971) 

103 

(117) 

11  (19) 

0.055 

(0.062) 

Waltham,  MA^ 

2/77-1/79 

960 

(971) 

114 

(117) 

17  (19) 

0.060 

(0.062) 

2 

Rosraan,  NC 

CY74 

1168 

(1149) 

127 

(119) 

52  (46) 

0.068 

(0.063) 

Austin,  TX 

2/78-1/79 

819 

(823) 

89 

(82) 

17  (41) 

0.047 

(0.044) 

Crawford  Hill, 

NJ3  CY68 

907 

(981) 

106 

(115) 

24  (25) 

0.056 

(0.061) 

Tampa ,  FL 

1/78-12/79 

1401 

(1254) 

107 

(107) 

82  (88) 

0.057 

(0.057) 

Slough,  ENG4 

6/68-6/72 

678 

(605) 

170 

(167) 

*  (15) 

0.092 

(0.089) 

1.  Boston  (Logan  Airport),  MA 

2.  Average  of  Greenville,  SC,  and  Asheville,  NC 

3.  Newark  (Newark  Airport),  NJ 

4.  Taken  from  Gatwlck  Airport  data  and  scaled  by  other  known  climatological  data 
to  Heathrow  Airport  (e.g.,  nominal  total  rainfall  at  Gatwick  is  1.13-times 
that  at  Heathrow) 

♦Unavailable  -  Use  nominal  value 
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Fig.  2.3.1.  Rain-rate  distributions  for  four  representative 
climate  areas. 
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times  were  dropped  to  0.1  second  or  less  it  is  expected  that  the  probability 
values  at  the  high  rates  will  increase  and  a  second  lognormal  distribution 
will  be  uncovered*. 

It  is  tempting  to  consider  a  crude  separation  into  two  types  of  rain, 
thunderstorm-created  and  non-thunderstorm-created.  If  the  long-term 

cumulative  statistics  for  rain-rates,  R,  are  defined  as  Pt(R>Rq)  for 
thunderstorms  and  Pnt(R>Rq)  for  non-thunderstorm  rain,  then  a  combined 
cumulative  statistic  might  be  written  as 

P(A  >  Aq)  =  [Dt/Dr]  Pt(R  >  Rq)  +  [1  -  Dt/Dr]  Pnt(R  >  Rq) 

(2.3.1) 

where  DT  is  the  number  of  thunderstorm  days  and  DR  is  the  number  of  days  of 
measurable  rain  as  defined  in  Section  2.2.  By  using  the  Portland,  OR,  and 
New  Orleans,  LA,  data,  two  lognormal  distributions  representing  Pt(R<Rq)  and 
Pnt(R>Rq)  are  generated  and  plotted  in  Fig.  2.3.2.  Their  fit  to  the  long-term 
measured  point  rain-rate  data  is  illustrated  in  Fig.  2.3.3. 

The  agreement  in  Fig.  2.3.3  provides  some  further  credibility  to  the 
empirical  determination  given  in  Section  2.2  and  indicates  that  the 

ratio  of  thunderstorm  days  to  days  of  measurable  rain,  DT/DR,  is  likely  to  be 
an  important  parameter  in  models  for  EHF  satellite-communications  systems 
behavior  with  rain.  It  should  be  noted,  however,  that  the  correlation  between 
point  rain  rate  and  attenuation  is  not  perfect. 

The  attenuation  from  rain  is  an  area  phenomena.  A  relation  between 
maximum  rainfall,  area  over  which  it  is  observed,  and  the  duration  was 
obtained  by  Fletcher  [1950] .  The  relation  becomes 

173 

Rt  *  12.7/T^  [1  +  -  ]  mm  (2.3.2) 

km 

ff 

In  fact,  the  maximum  point  rain  rates  are  bounded  by  the  reciprocal  of  the 
square-root  of  their  duration  [Fletcher,  1950]  and  a  decrease  in  averaging 
time  by  over  two  orders-of-magnitude  should  provide  an  additional  order-of- 
magnitude  in  observed  long-term  rain-rates. 
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PROBABILITY  RATE  EXCEEDS  INDICATED  VALUE 


PROBABILITY  RATE  EXCEEDS  INf  iCATED  VALUE 


Fig.  2.3.3.  Reconstructed  rain-rate  distributions  using  data  from 
Fig.  2.3.2.  Measured  data  points  are  given  for  comparison. 
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where 


ftj,  is  the  total  rainfall  in  am 
Tq  is  the  rainfall  duration  in  hours 
Dkm  is  the  diameter  in  km  of  the  area  over  which  the  rain  amount,  Rt, 
falls  in  time  TR 

Note  that  a  point  rainfall  maximum  envelope  becomes  (0^=0) 

, POINT  "  364/^  "  <2-3'3> 

and  a  wide-area  rainfall  maximum  envelope  becomes  (Djtm=*co) 

“t.uiob  -  12-7^  -  <2-3-4> 

If  the  "effective”  rain  rate  becomes  R-Rj/Tq  mm/hr,  then  it  should  be  expected 
that  the  most  intense  rainfalls  should  last  for  the  shortest  times  and  extend 
over  the  smallest  areas.  This  factor  has  obvious  application  in  diversity 
siting  as  well  as  "waiting-time"  determinations. 

2.4  Elevation  Angle  Dependence 

Because  the  path  lengths  in  the  atmosphere  are  different,  it  would  be 
expected  that  two  different  elevation  angles  on  the  same  azimuth  examined 
simultaneously  during  a  rainstorm  would  exhibit  different  values  of 
attenuation.  At  low  rain  rates,  where  rainfall  is  generally  widespread  and 
relatively  homogeneous  both  in  azimuth  and  elevation  (up  to  the  0°-C 
isotherm),  the  isochronal  attenuations  would  be  approximately  related  by  the 
classical  csc8-relatlonship*.  However,  the  more  intense  rain,  particularly 
that  originating  in  thunderstorms,  is  by  no  means  homogeneous,  either  in 
azimuth  or  in  elevation. 


*0r  modifications  thereto  caused  by  a  spherical  earth. 
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Even  simple  geometries  like  the  cylindrical  rain  cell  illustrated  In 
Fig.  2.4.1  do  not  yield  simple  results.  For  example,  assume  the  rain  extends 
uniformly  to  the  ground  from  the  melting  layer  height,  hg  (0°-C  Isotherm). 
Even  If  the  edge  of  the  rain  cell  Includes  the  terminal  (R-d/2),  the  csc0 
behavior  will  only  be  apparent  if  d»hg.  If  d<hg,  as  can  be  the  case,  for 
much  of  the  useful  elevation-angle  range  the  attenuation  due  to  that  rain  cell 
increases  with  increasing  0.  Other  more  complex  geometries  will  give  more 
complex  elevation-angle  dependencies. 

The  yearly  behavior  of  hg  at  two  different  locations  (obtained  from 
National  Summaries  of  Local  Climatological  Data)  is  shown  in  Fig.  2.4.2.  Note 
that  melting- layer  heights  in  excess  of  4  km  are  common  for  those  months  of 
the  year  in  which  thunderstorms  are  prevalent.  Thunderstorm  rain-cell  sizes 
measured  by  a  weather  research  radar  in  Goodland,  Kansas  [Crane,  et  al.,  1977; 
Crane,  1980]  are  converted  by  the  relationship  for  thunderstorm  rain  rate,  R, 
and  radar  reflectivity,  Z  [Rogers,  1976]  as 

Z  -  200R1*6  .  (2.4.1) 

The  data  are  sparse,  but  an  approximate  lognormal*  behavior  is  obtained  and 
Illustrated  in  Fig.  2.4.3.  Note  that  as  rain  rate  (and  potential  attenuation) 
increases,  the  median  cell  diameter  drops  from  about  2.5  km  at  low  rain  rates 
(~0.5  mm/hr)  to  about  1.6  kra  at  moderate  to  high  rain  rates  (~40  mm/hr). 
Therefore,  most  of  the  time,  thunderstorm  rain  cells  will  have  effective 
diameters  much  smaller  than  their  heights  and  the  csc-0  behavior  for 
attenuation  will  not  be  observed. 

In  many  cases,  not  only  are  the  rain-cell  structures  taller  than  their 
effective  diameters,  but  they  are  tilted  and  distorted  by  the  natural 
progression  velocity  of  the  storm  [Goldhirsh,  1976].  Furthermore,  there  is  no 
guarantee  that  the  rain  rates  within  a  given  cell  will  be  homogeneous  with 
altitude;  such  a  lack  of  vertical  uniformity  is  illustrated  by  Fig.  2.4.4, 

*Similar  lognormal  size  behavior  has  been  observed  by  Strickland  [1974). 


Fig.  2.4.1.  Illustration  of  a  simple  cylindrical  rain-cell  geometry 
with  melting-layer  height  h  ,cell  diameter,  d,  distance  from  center 
of  cell  to  terminal  of  R  an3  satellite  elevation  angle,  0. 


(km) 


Fig.  2.4.2.  Average  zero-degree  (C)  isotherm  heights  as  a  function  of 
time  over  a  year  for  two  representative  locations. 
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j  Fig.  2. A. 3.  Cumulative  probability  distribution  functions  for 

rain-cell  diameters  with  rain  rate  as  a  parameter. 
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taken  from  the  Handbook  of  Geophysics  [1961]  and  based  on  data  obtained  by 
Atlas  [1948].  A  terminal  located  at  the  zero-kilometer  reference  point  and 
looking  east  through  the  storm  would  experience  a  peak  path  rain  rate  of 
288  mm/hour  at  9=30°  dropping  to  70  mm/hour  at  6=7°.  It  also  appears  that  the 
cell  structure  is  such  that  the  total  path  integrated  rain  rate  at  9=30°  would 
be  greater  than  that  observed  at  9=7°  by  a  factor  even  larger  than  the  ratio 
of  peak  path  rain  rates. 

Finally,  the  atmospheric  portion  of  the  earth-space  path  could  well 
encompass  several  different  rain  cells  at  the  same  time.  A  transverse  cross- 
section  through  a  thunderstorm  [Crane,  1971,  and  Crane  and  Blood,  1979]  is 
given  in  Fig.  2.4.5.  At  azimuths  of  279.5°  and  282.5°  more  than  one  rain  cell 
would  be  present  in  the  earth-space  path  simultaneously.  As  the  storm  moves, 
the  complex  structure  of  rain  cells  will  form,  change  and  reform  while  the 
earth-space  path  "samples"  the  non-stationary  behavior. 

This  azimuthal  complexity,  coupled  with  the  non-uniform  distribution  with 
elevation  and  the  small  dlameter-to-height  ratios  of  typical  intense  rain 
cells,  makes  the  determination  of  any  elevation-angle  dependence  of  path 
attenuation  very  difficult.  Consequently,  statistics  gathered  on  path 
attenuation  over  long  periods  of  a  year  or  more  ought  to  have  a  complex 
mixture  of  all  forms  of  path  and  elevation-angle  dependences.  It  seems 
likely,  therefore,  that  such  long-term  statistics  should  show  little  or  no 
elevation-angle  dependence,  and  models  for  system  availability  taken  from  such 
statistical  data  should  not  show  elevation-angle  dependence*. 

2.5  Some  Observations  on  the  Time  Structure  of  Rain-Induced  Fading 

Although  the  measures  used  in  this  report  are  based  on  1-minute  period 
rain  rates,  much  more  rapid  fluctuations  in  point  rain  rates  exist.  Freeny 
and  Gabe  [1969]  found  rapid  second-to-second  fluctuations  in  fractional-second 
rain  rates  of  up  to  3  to  1  with  8-inch  diameter  "instantaneous"  rain  gauges. 

"kit 

Since  rain  cells  move  at  peak  velocities  of  ~14  m/sec  [Englebrecht ,  1979] 


Much  of  the  data  used  in  this  report  was  taken  from  radiometric  sun-trackers 
and,  therefore,  represents  a  mixture  of  data  from  a  continuum  of  elevation 
gngles.  Any  elevation-angle  dependence  would  not  be  retrievable. 

Average  velocities  are  much  lower,  l  to  4  m/sec  [Cra'.e,  1980]. 


RANGE  (km) 


TRANSVERSE  DISTANCE  AT  150km  (km) 


Fig.  2.4.5.  Radar-derived  range-azimuth  cut  at  constant  elevation 
angle  (1.4°)  through  a  thunderstorm. 
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and  fine  structure  of  the  order  of  1  m  [Lin,  1975]  in  the  cells  can  be  moved 
by  winds  over  100  kms  (~28  m/sec),  such  behaviors  are  to  be  expected. 
However,  they  are  not  necessarily  representative  of  the  path-integrated  rain- 
rate  behavior  responsible  for  the  total  path  attenuation  at  any  given  time. 

On  short-haul  terrestrial  paths  of  length  L  the  effective  volume  of  the 
path  over  which  the  rainfall  will  be  averaged  is  the  Fresnel  ellipsoid  of 
major  axis  L  and  minor  axes  /XL  [Bodtmann  and  Ruthroff,  1974].  For  an 
earth/satellite  path,  however,  the  effective  averaging  volume  is  encompassed 
by  a  paraboloid  directed  along  the  path  as  illustrated  in  Fig.  2.5.1.  The 
height,  hQ,  of  the  melting  layer  and  the  elevation  angle,  9,  determine  the 
effective  rain-containing  path  length,  L.  The  radius  of  the  Fresnel  zone  is 
/XL  and  the  volume  is  irXL^/2. 

If  hQ  is  set  as  the  typical  maximum  height  of  the  melting  layer  of  5  km 
(perhaps  higher  inside  large  thunderstorms)  and  9=30°,  the  path  length 
L=10  km.  For  a  wavelength  of  0.01  m  (30  GHz),  the  Fresnel  radius  at  the  far 
end  of  the  path  is  only  10  ra  and  the  volume  is  1.57x10  m  .  For  light  rain  of 

O 

0.25  mm/hr  there  are  only  about  20  drops/m  [Lin,  1975]  and  the  number  of 
drops  in  the  Fresnel  zone  would  only  3.14x10^. 

The  averaging  effect  of  the  Fresnel  zone,  both  in  width  and  in  length, 
restricts  the  time  rate  of  change  of  attenuation.  A  "square"  rain  cell  moving 
at  14  m/sec  normal  to  the  path  would  fill  the  path  volume,  if  it  were  located 
near  the  melting  height  in  the  previous  example,  in  about  0.7  second.  Fine 
structure  of  the  order  of  1  m  would  transit  the  path  in  less  than  1/4  of  a 
second  if  being  moved  by  100-kph  winds  normal  to  the  path.  The  fine  structure 
would,  however,  contribute  far  less  to  the  attenuation  than  would  the  cell 
itself.  Consequently,  the  attenuation  could  be  expected  to  change,  in  the 
worst  case,  several  dB  in  a  few  seconds  (1  to  3)  with  a  low-level  "noise"  on 
the  attenuation  with  spectral  components  constrained  to  the  region  of  a  few 
Hertz  (1  to  5)  or  less  . 

•Jg 

Their  geometry  would  increase  the  attenuation  "rise  time"  for  cylindrical 
cells . 
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Fig.  2.5.1.  First  Fresnel  zone  for  a  terminal-satellite  path  of 
elevation  angle,  0 ,  to  a  melting  layer  height,  h  ,  with  consequent 
effective  path  length,  L.  ° 
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Fig.  2. 6. 2.1.  Attenuation  data  for  Crawford  Hill,  NJ. 
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Fig.  2. 6. 2. 2.  Fade-duration  distributions  for  Crawford  Hill,  NJ. 
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Fig.  2. 6. 2. 3.  Margin  Increase  for  Crawford  Hill,  NJ. 
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Fig.  2. 6. 3. 2.  Fade-duration  distributions  for  Slough,  ENG,  at  19GHz. 
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From  these  numerical  examples,  it  can  be  seen  that  the  onset  of  rain- 
induced  fades  can  be  rapid,  but  not  instantaneous.  A  few  seconds  should  be 
available  to  sense  the  onset  of  the  fade  and  to  take  corrective  action. 
Although  there  will  be  fine  structure  on  the  attenuation  with  significant 
spectral  components  below  a  few  Hertz,  such  fine-structure  amplitude 
modulation  should  only  affect  those  systems  where  significant  carrier-to-nolse 
ratios  below  a  few  Hertz  from  the  carrier  are  required.  It  is  interesting  to 
note  that  these  "rapid"  fine-structure  modulations  would  still  only  amount  to 
a  spectral  spread  of  a  few  parts  in  10^  at  30  GHz  and  would  probably  be 
overwhelmed  by  the  short-term  noise  of  the  frequency  standards  used  in  the 
terminal  and/or  satellite. 

Finally,  some  of  the  data  used  in  this  report  was  obtained  from  sun-track 
experiments.  The  sun  is  not  a  point  source  but  instead  occupies  an  angular 
extent  of  about  8.7  milliradlans.  At  30  GHz  the  10-km  Fresnel  zone  occupies 
an  apparent  angle  of  2  milliradlans.  Therefore,  the  high  frequency  time 
structure  of  attenuation  observed  by  a  sun  tracker  might  be  expected  to  be 
somewhat  different  than  that  observed  on  a  satellite  downlink.  This 
difference  occurs  for  much  the  same  reason  that  stars  visibly  "twinkle"  under 
natural  atmospheric  turbulence  but  planets  do  not.  Some  additional  averaging 
of  path  attenuation  fluctuations  is  present,  if  the  source  has  significant 
angular  extent. 

2.6  Model  Inputs 

2.6.1  Data  Ordering 

Each  group  of  data  from  a  given  site  and  given  frequency  is  reduced,  if 
necessary,  to  the  distribution  function  of  attenuation  when  it  is  raining, 
with  the  aid  of  Pg^IN  38  detertnined  in  Section  2.3.  From  the  P(t>tq)  and  the 
P(A>AQ|given  rain)  the  relation  in  (2.1.7)  will  give  a  value  of  M2  for  each  M^ 
and  storage  time,  T-tq.  All  of  the  data  points  are  gathered  together  on  one 
plot  of  Mj/M^  versus  T. 

The  points  on  the  individual  plots  of  M2/M^  versus  T  are  fitted  to  the 
empirical  relation 
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(2. 6. 1.1) 


m2/m1  -  (1  +  |-)0'5 

o 

with  the  value  of  TQ  chosen  to  give  the  best  mean-square-error  fit  to  the 
data.  It  Is  Interesting  to  note  that,  with  few  exceptions,  the  data  points  on 
the  M2/M^  versus  T  plots  show  little  or  no  variation  with  or  operating 
frequency.  There  does  appear  to  be  a  variation  In  TQ  from  site  to  site, 
however.  The  reasons  for  such  variation  will  be  examined  In  Section  2.7. 

2.6.2  Crawford-Hlll,  NJ,  Data 

Sun-tracker  experiments  were  carried  out  by  Bell  Telephone  Laboratories  at 
Crawford  Hill,  NJ,  from  12/67  to  2/69  [Wilson,  1969,  and  Lin,  1973].  The 
16-GHz  and  30-GHz  distribution  functions  for  P(A>AQlgiven  rain)  are  given  In 
Fig.  2. 6. 2.1  and  the  distribution  functions  for  duration  of  fade,  P(t>tq),  for 

9- dB  and  21-dB  thresholds  at  30  GHz  and  a  9-dB  threshold  at  16  GHz  are  given 
in  Fig.  2. 6. 2. 2. 

It  is  Interesting  to  note  that  the  data  presented  In  Figs.  2.6.2. 1  and 
2. 6. 2. 2,  plotted  on  "log-normal-straightening"  paper,  fit  straight  lines 
rather  well.  The  only  exception  appears  at  the  long  durations  at  the  lowest 
frequency.  Such  deviations  are  to  be  expected  In  a  finite  data  set  [Lin, 
1976b] ,  where  the  tails  of  the  distributions  would  not  be  expected  to  be 
precisely  lognormal. 

The  generated  plot  of  f^/M^  versus  T  is  given  in  Fig.  2. 6. 2. 3.  A  value 
for  Tq  of  2.6±0.7  minutes  fits  the  data  best  and  the  resulting  curve  of 
(2.6. 1.1)  is  given  on  the  same  figure. 

2.6.3  Slough,  GNG,  Data 

Sun-tracker  experiments  were  carried  out  by  the  GPO  at  Slough,  ENG,  at 
19  GHz  from  9/68  to  7/72  and  at  37  GHz  for  the  rather  brief  period  of  9/71  to 
6/72  [Davis,  1971,  and  Davis,  1973].  The  resulting  curves  of  PCAMglgiven 
rain)  are  given  In  Fig.  2. 6. 3.1.  Fade-duration  distributions  for  5-dB  and 

10- dB  thresholds  at  19  GHz  are  given  in  Fig.  2. 6. 3. 2  and  at  38  GHz  in 
Fig.  2. 6. 3. 3. 

Again,  a  good  fit  to  lognormal  distributions,  except  at  the  extreme  tails, 
is  noted.  The  data  points  in  Fig.  2. 6. 3. 4  for  h^/M^  versus  T  are  surprisingly 
tight  and  well-fitted  by  TQ“4.15±0.97  minutes. 
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2. 6. A  Rosraan ,  NC,  Data 

Satellite-link  experiments  were  carried  out  at  30  GHz  from  the  NASA  site 
at  Rosman,  NC  [Ippolito,  1975,  and  Crane,  1980].  The  period  covered  only  a 
few  of  the  months  (summer)  of  197A  so  the  data  are  not  expected  to  be  a  good 
representation  of  a  year-long  sample.  However,  P(A>AQ|given  rain)  is  plotted 
in  Fig.  2.6.4. 1  and  the  fade-duration  distributions  for  5-dB,  10-dB  and  15-dB 
thresholds  are  given  in  Fig.  2. 6. 4. 2. 

The  resulting  M2/M^  behavior  versus  T  is  given  in  Fig.  2. 6. 4. 3.  As 

expected,  the  scatter  is  large  and  the  rather  loose  fit  is  also  implied  by  a 
value  of  T0“l .42±1 • 17  minutes. 

2.6.5  Waltham,  MA,  Data 

Satellite-link  experiments  using  CTS  were  carried  out  at  11.7  GHz  at 
Waltham,  MA,  over  the  two  periods  2/77  to  1/78  and  2/78  to  1/79  [Nackoney, 

1979,  and  Tang,  1979].  The  resulting  P(A>AQlgiven  rain)  is  given  in 

Fig.  2.6. 5.1.  Fade-duration  distributions  for  3-dB,  6-dB  and  10-dB  thresholds 
are  given  in  Fig.  2. 6. 5. 2.  The  M2/M^  versus  T  plot  of  Fig.  2. 6. 5. 3  yields  a 
value  of  T»2.08±1.53  minutes. 

Note  that  the  two  years  of  data  in  Fig.  2.6. 5.1,  if  separated,  show 
distinct  abnormalities  at  the  low-probability/high-attenuation  end  of  the 
curves.  The  period  2/78  to  1/79  deviates  significantly  from  lognormal  in  such 
a  fashion  that  high  attenuation  events  would  be  expected  to  be  more  rare  for 
that  period  than  for  2/77  to  1/78.  In  fact,  as  shown  in  Table  2.3.2  in 

Section  2.3,  the  number  of  thunderstorm  days  for  the  period  2/78  to  1/79  is 
11,  one-half  of  the  value  of  22  for  2/77  to  1/78  and  about  one-half  of  the 
average  value  of  19.  Therefore,  the  lack  of  a  normal  number  of  high- 
attenuation  events  for  the  period  2/78  to  1/79  is  expected;  however,  the 
average  of  19  thunderstorm  days  for  the  two-year  period  is  close  to  nominal 
and  the  composite  data  is  consequently  close  to  the  expected  lognormal 
behavior. 

2.6.6  Austin,  TX,  Data 

The  11.7-GHz  signal  from  the  CTS  satellite  was  also  used  to  gather 
attenuation  data  at  Austin,  TX,  for  the  period  2/78  to  1/79.  The  plot  of 
PCA^A^jIgiven  rain)  is  illustrated  in  Fig.  2.6.6. 1  and  3-dB,  6-dB  and  10-dB 
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2.6.4. 1.  Attenuation  data  for  Rosman,  NC. 
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Fig.  2. 6. 5. 3.  Margin  increase  for  Waltham,  MA. 
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threshold  fade-duration  distributions  are  given  in  Fig.  2. 6. 6. 2.  A  value  of 
Tq*1.85±0.63  minutes  is  indicated  by  the  versus  T  behavior  in 

Fig •  2.6. 6. 3. 

As  in  the  second  year  of  the  Waltham,  MA,  data,  there  is  a  significant 
lack  of  high-attenuation  events  in  the  Austin,  TX,  data.  From  Table  2.3.2  the 
number  of  thunderstorm  days  for  the  observing  period  was  only  17,  less  than 
one-half  the  nominal  value  of  41.  Therefore,  the  period  was  abnormal  and  a 
data  take  over  several  years  would  be  expected  to  have  a  better  distribution 
in  the  tails. 

2.6.7  Tampa,  FL,  Data 

Data  were  taken  over  the  period  1/78  to  12/79  using  the  COMSTAR  19-GHz 

beacon  at  3  sites  in  the  Tampa,  FL,  area  [Tang,  1979,  and  Tang,  1980].  The 
complete  attenuation  data  for  the  University  of  Southern  Florida  (USF)  site 
were  supplied  by  Tang  [1980]  as  a  personal  communication  and  are  plotted  in 
Fig.  2. 6. 7.1.  The  fade-duration  data  are  given  for  6-dB  and  12-dB  thresholds 
in  F ig .  2. 6. 7.2. 

The  fade-duration  data  exhibit  a  strange  behavior  at  both  the  "low- 
probability"  tails  and  the  "high-probability"  tails.  The  deviations  at  the 
low-probability  extreme  tails  can  be  explained  on  the  basis  of  a  limited 
sample  of  data.  The  high-probability  anomaly  is  a  different  matter.  It  is 

not  clear  if  it  results  from  an  instrumentation  or  computational  "glitch"  or 

if  it  is  real*. 

Because  of  this  anomaly  two  different  I^/M^  versus  T  data  sets  are 

presented  in  Fig.  2. 6. 7. 3.  One,  termed  "actual",  follows  the  data  points  in 
Fig.  2. 6. 7. 2.  The  second,  termed  "projected"  follows  the  projected  lognormal 
behavior  in  Fig.  2. 6. 7. 2.  There  are,  therefore,  two  values  of  T  ,  a 
"projected"  value  of  To“0.43±0.13  minute  and  an  "actual"  value  of 
TQ*0.29t0.05  minute. 


*If  there  were  a  slow-response  problem  in  the  amplitude-measuring  equipment, 
the  shorter  fades  might  well  be  "filled  in”  and  missed. 
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Fig.  2. 6. 6. 2.  Fade-duration  distributions  for  Austin,  TX. 
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Each  of  the  data  sets  from  different  locations  was  analyzed  in  Section  2.6 
and  a  storage-time  relationship  of  the  form 

M2/M1  =  [1  +  |~]°'5  (2.7.1) 
o 

was  derived  for  each  location.  The  values  of  TQ,  however,  varied  from  a  low 
of  0.43  minute  in  Tampa,  FL,  to  a  high  of  4.15  minutes  in  Slough,  ENG.  It 
remains  to  determine  if  certain  climatological  factors  can  be  identified  as 
the  causes  of  this  nearly  ten-to-one  change  in  the  parameter  TQ. 

The  characteristics  of  rainfall  from  thunderstorms  are  dominated  by  short- 
duration  high-rate  events,  whereas  the  more  widespread  rainfall  is  more  likely 
to  have  long-duration  moderate-rate  events.  Therefore,  it  should  be  expected 
that  the  storage-time  requirements  for  Tampa,  FL,  where  much  of  the  rainfall 
is  from  thunderstorms,  would  be  less  than  those  for  Slough,  ENG,  where  most  of 
the  rainfall  is  not  thunderstorm-derived.  As  noted  in  Section  2.3,  a  useful 
index  is  the  ratio  of  thunderstorm  days  to  days  of  measurable  precipitation, 

dt/dr. 

Values  of  Dt/Dr  for  the  periods  of  time  and  locations  are  taken  from 
Section  2.3  and  paired  with  the  TQ  values  from  Section  2.6  in  Table  2.7.1. 
There  are,  of  course,  many  different  function  behaviors  for  T0(dt/°r)  that  can 
be  proposed.  However,  by  plotting  logTQ  versus  Dt/Dr  in  Fig.  2.7.1,  a  simple 
straight-line  model  can  be  deduced. 

The  model  agreement  is  improved  somewhat  by  realizing  that  the  sun-tracker 
experiments  at  Crawford  Hill,  NJ,  and  Slough,  ENG,  typically  do  not  operate 
when  the  sun  is  less  than  20°  above  the  horizon.  Consequently,  they  will  miss 
the  early-evening  and  night-time  population  of  thunderstorms.  It  is  somewhat 
arbitrarily  assumed  that  roughly  one-third  of  the  thunderstorms  will  be  missed 
by  such  operation  and  the  effective  value  of  Dt/Dr  will  be  reduced  to  two- 
thirds  of  its  original  value.  Such  a  reduction  is  shown  by  the  "dashed" 
circles  in  Fig.  2.7.1. 


TABLE  2.7.1 

VALUES  OF  TQ  AND  THUNDERSTORM  INDEX,  Dt/Dr  FOR  THE  VARIOUS  LOCATIONS 


LOCATION 

dt/dr 

T  ,  minutes 

Crawford  Hill,  NJ 

0.226 

2.6  ±  0.7 

Slough,  ENG 

0.088 

4.15  ±  0.97 

Rosman,  NC 

0.409 

1.42  +  1.17 

Waltham,  MA 

0.149 

2.08  ±  1.53 

Austin,  TX 

0.191 

1.85  ±  0.63 

Tampa ,  FL 

0.766 

(actual)  0.29  ±  0.13 
( projected)  0.43  ±  0.05 
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From  Fig.  2.7.1,  the  model  value  for  TQ  is 

-3.7D_/D 
T  R 

T  3  5e  minutes 

o 


(2.7.2) 


giving  a  storage-time  relation  of 


r  3.7D  /D  -10.5 

m2/m1  3  L1  +  °-2T  e  j 


(2.7.3) 


Since  1>Dt/Dr>0,  the  relationship  in  (2.7.3)  is  bounded  so  that  climate  areas 
like  the  Northwest  coast  of  the  US  with  Dt/Dr  0  and  other  areas  like 
equatorial  Africa  and  equatorial  South  America  with  Dt/Dr  1  are  both 
represented.  There  is,  therefore,  some  reason  to  believe  that,  within  the 
limits  of  the  data,  the  empirical  storage-time  model  for  (2.7.3)  is 
geographically  universal. 

Finally,  some  measure  of  the  scatter  in  the  data  can  be  determined  by 
plotting  M2/Mj  versus  the  normalized  storage  time,  T/Tq,  for  all  of  the 
locations  on  one  graph.  Such  a  plot  is  given  in  Fig.  2.7.2.  Note  that  the 
four  anomalous  points  on  the  low  side  of  the  composite  come  from  the  15-dB 
threshold  values  for  Rosman,  NC,  and  represent  an  inadequate  data  base.  There 
is  also  some  reason  to  believe  that  threshold  values  of  3  dB  or  less  do  not 
yield  good  results.  The  3-dB  threshold  points  from  various  locations 

represent  many  of  the  widely  scattered  anomalous  points  on  the  high  side  of 
the  composite  curve.  With  such  "bad”  data  excluded,  over  a  range  of 
6>M2/M^>1,  the  required  value  of  T/TQ  determined  directly  from  the  data  is 
less  than  a  factor  of  1.5  different  than  the  predicted  (model)  value. 

2.8  An  Example 

In  order  to  illustrate  the  use  of  the  waiting-time  model,  consider  a 
satellite  system  with  a  terminal  in  Crawford  Hill,  NJ  at  an  uplink  frequency 
of  30  GHz*.  By  combining  the  nominal  value  of  PraIN”^’®^  (from  Section  2.3) 


A  44-GHz  example  might  be  prepared  by  using  the  frequency  dependence 
illustrated  by  the  NJ  and  ENG  data  and  extrapolating  30-GHz  data  to  44  GHz. 


m2/m, 


for  the  location  with  the  lognormal  distribution  for  the  attenuation  (from 
Section  2.6),  a  yearly  availability  curve  can  be  generated.  Note  that  the 
resulting  availability  function  as  given  in  Fig.  2.8.1  is  not  lognormal*. 

For  the  nominal  values  of  D^*24  and  D^=106,  from  Section  2.3,  the  storage¬ 
time  relationship  from  (2.7.3)  becomes 

M2/M1  -  [1  +  0.462T]0,5  (2.8.1) 

for  T-minutes  of  storage.  The  resulting  improvement  in  availability  with 
T*l,3,10,30  and  100  minutes  is  illustrated  on  Fig.  2.8.1.  The  curves  are  not 
continued  much  below  3  dB  because  it  is  expected  that  the  model  becomes 
progressively  more  Inaccurate  as  the  initial  margin  becomes  less  than  3  dB. 

Note  that  slightly  less  than  20  minutes  of  storage  will  produce  an  order- 
of-magnitude  increase  in  availability  over  a  wide  range  of  initial  margins. 
At  25  dB  of  margin,  for  example,  the  availability  is  increased  from  0.999  to 
0.9999  by  supplying  the  system  with  20-minute  data-storage  capability.  The 
same  increase  in  availability  without  storage  would  require  increasing  the 
margin  from  25  to  about  80  dB.  At  lower  initial  margins  somewhat  more  storage 
is  required  to  produce  the  or der-of -magnitude  increase  (e.g.,  at  6  dB  for 
over  30  minutes  of  storage  is  necessary). 

Although  it  would  seem  that  an  availability  producing  an  average  outage  of 
1  hour  per  year  (~0.9999)  would  reach  an  availability  of  1.0  with  1  hour  of 
storage,  such  is  not  so.  The  1  hour  figure  is  for  an  average  year  and  there 
is  still  a  quite  finite  probability  that  the  outage  could  be  well  in  excess  of 
1  hour  if  it  was  only  likely  to  occur  once  every  several  years.  Therefore, 
the  100-minute  curve  in  Fig.  2.8.1  does  not  produce  a  unity  availability. 


A 

Obviously,  the  availability  must  be  greater  than  1-P^^  re8ardless  of  how 
hard  (or  long)  it  may  rain. 
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III.  CONCLUDING  REMARKS  -  FUTURE  WORK 

Utilizing  the  rather  limited  attenuation  data  available,  a  model  for  the 
effect  of  storage  time  on  rain-induced  EHF  availability  can  be  developed.  The 
model  appears  to  be  "universal"  and  seems  applicable  to  virtually  any 
climatologic  region.  The  parameters  which  generate  the  model  are  readily 
available  from  the  standard  climatological  data  for  any  location  on  the  globe 
where  such  data  is  regularly  gathered. 

The  improvement  in  effective  margin  can  be  written  as 

3.7  D  /D  0.5 

M2/M1  =*  1  +  0.2  T  e  (3.1) 


where 

T  is  the  data  storage  time  in  minutes, 

D,j,  is  the  number  of  thunderstorm  days/year, 

Dr  is  the  number  of  days/year  of  rain  (defined  as  the  number  of 
days/year  of  precipitation  over  0.25  mm  minus  the  number  of 
days/year  of  snowfall  greater  than  25. A  mm) 
is  the  link  margin  in  dB  with  no  storage,  and 
M2  is  the  equivalent  link  margin  in  dB  with  storage  of  T  minutes 
available 


It  seems  reasonable  to  expect  that  additional  data  would  result  in  some  change 
in  the  parameter  characterizing  (3.1),  but  the  general  behavior  should  remain 
unchanged . 

The  determination  of  long-term  attenuation  distributions  and  fade-duration 
distributions  is  a  difficult  and  time-consuming  task.  Much  of  the  data 
gathering  and  manipulation  could,  however,  be  automated  and  accumulated 
virtually  in  real  time,  if  appropriate  hardware  were  developed  and  distributed 
to  many  of  the  satellite  ground  stations.  More  EHF  beacons  should  be  packaged 
onto  available  platforms.  If  such  a  program  were  implemented,  it  would  be 
necessary  to  choose  an  ensemble  of  sites  which  provide  a  good  sample  of  the 
various  climatologic  regions,  in  particular  those  regions  with  a  very  low 
incidence  of  thunderstorm  activity. 
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